Surface bound states and spin currents in non-centrosymmetric superconductors 
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We investigate the influence of spin-orbit coupling in a non-centrosymmetric superconductor on its 
ground state properties near a surface. We determine the spectrum of Andreev bound states due to 
surface- induced mixing of bands with opposite spin helicities for a Rashba-type spin-orbit coupling. 
We find a qualitative change of the Andreev spectrum when we account for the suppression of the 
order parameter near the surface, leading to clear signatures in the surface density of states. We 
also compute the spin current at the surface, which has spin polarization normal to that of the bulk 
current. The magnitude of the current at the surface is enhanced in the normal state compared to 
the bulk, and even further enlarged in the superconducting phase. The particle and hole coherence 
amplitudes show Faraday-like rotations of the spin along quasiparticle trajectories. 

PACS numbers: 72.25.-b, 74.45. +c, 73.20.At 
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The role of chirality and spin-orbit coupling in ma- 
terials and nanostructures is a very active and promis- 
ing subject in the fields of spintronics, superconductiv- 
ity and magnetism [l|, l2, Isj . The unusual properties of 
non-centrosymmetric (NCS) materials originate from the 
crystal structure that lacks a center of inversion, allow- 
ing for pronounced spin-orbit (SO) coupling that is odd 
in the electron momentum, and leading to a chiral ground 
state. The resulting two-band nature of NCS metals 
leads to effects reminiscent of semiconductor physics, 
such as birefringence and spin polarization of the electron 
wavepacket [j]. Especially promising is the presence of 
charge-neutral spin currents in the ground state [5|, l6| . 

Since understanding of interface physics is one of the 
foundations for all potential applications, it is of pivotal 
interest to investigate how the physical properties of NCS 
materials are modified near surfaces. The key observa- 
tion is that scattering events off interfaces in materials 
with strong spin-orbit effects are typically spin-active. 
Spins dominate the surface physics, and any successful 
theoretical treatment must take this into account. 

The recently discovered class of NCS superconductors 
[71, [sl, Isl, uM combines the strong SO coupling that gov- 
erns the metallic bands with a non-trivial, chiral, spin 
structure of the superconducting (SC) order parameter 
12l.ll3|. As a result, one may expect that spin trans- 
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port in the SC phase exhibits novel features compared to 
superconductors with negligible SO interaction. These 
features are expected to be especially prominent near 
surfaces and interfaces, where the physics is controlled 
by the Andreev bound states (ABS), built as a result of 
particle-hole coherent scattering. ABS are crucially im- 
portant in unconventional superconductors [lj|, where 
the phase variation of the order parameter (OP) on the 
Fermi surface [l5| and the pairbreaking near interfaces 
may lead to a midgap peak in the density of states (DOS) 
at the surface. The ABS states control thermodynamic 



properties and stability of the surfa ce p hases 16|,ll7|, and 
govern transport across interfaces 18|, ll9|, |2C | . 

In this Letter we study the Andreev states and spin 
currents at the surface of a NCS superconductor. We 
show that (i) the ABS spectrum is qualitatively modified 
by the self consistent suppression of the order parameter; 
(ii) the spin current is (a) strongly enhanced near the 
surface in the normal state; (b) further enhanced in the 
SC phase. We develop a detailed theory of these effects. 

For a non-centrosymmetric material it is convenient 
to perform a canonical transformation from a spin basis 
(with fermion annihilation operators Ck^j. for spin /i =| , |) 
to the so-called helicity basis (&ks with hclicity s = ±), 
that diagonalizes the kinetic part of the Hamiltonian, 



Ukzn = Yl '^kM(^k + agkO-)Mi'Cki/ = Yl ^ksbisKs 



(1) 
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Here, ^k is the band dispersion relative to the chemical 
potential in the absence of SO interaction, a is the SO 
coupling strength, cr is the vector of Pauli matrices, and 
gk is a normalized (see below) SO vector llj, [l^ that is 
odd in momentum, g_k = — gk, see Fig. [TJ The helicity 
band dispersion is £k± = Ck ± Q;|gk|- SO interaction fixes 
the orientation of the quasiparticle spin with respect to 
its momentum in each hclicity band. 

The Hamiltonian ([T]) is time reversal invariant but lifts 
the spin degeneracy. The transformation from spin to 
hehcity basis, C/k, is defined by U^iSk^^Wl = Igkli^s, and 
determined by the direction of the g-vcctor in k space. 
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gk X z 

|gk X z| 



(2) 



where z is the unit vector in z-direction, and 6^ is the 
polar angle between gk and z |ia |. 

To describe superconductivity we use the Nambu- 
Gor'kov formalism modified for a helical basis. Wc define 
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FIG. 1: (Color online) A map of the spin-orbit vector in mo- 
mentum space for the Rashba form gk = k x z. On reflection 
the spin orbit vector gk niay change, e.g. from A ^r A' , or 
not, B — » B' . The scattering geometry is shown on the left. 



the helical counterpart, B^ — [b]^^, &](._, ^k+i ^k-)' ^^ t^*^ 
Nambu spinor, C^ = (4^ 4i' ^kt' Ck^), by 



velocity renormalization of order a/Ej ^ 1 is neglected. 
We choose a separable pairing interaction consistent with 
the form of the gap, and determine A self consistently 
with g. In NCS superconductors the OP is a mixture of 
spin singlet (Ag) and triplet (At) components |13l. l22j. 
Assuming that the triplet component aligns with gk, in 
real gauge it is given by, 

A = Kk/)[A,(R)i + At(R)z)so(k/)] {ia2)h , (5) 

where the basis function ^(k/) transforms according to 
one of the irreducible representations of the crystal point 
group, and (y^(k/)) = 1. With the gap functions in the 
helicity bands, A± = A^ ± At|gk^|, the order parameter 
is A = {A+(T+ - A_o-"}3^, where a^ = (cri ±i(j2)/2. 

We parameterize the Green's function by the coherence 
functions for particles and holes, 7 and 7 (2x2 spin matri- 
ces), which allow a very intuitive ph ysic al interpretation 
of the Andreev scattering processes [23| , 
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(6) 



Sk — f/kCk, U — 



U 

u* 



(3) 



and construct 4x4 retarded Green's functions in helicity 
basis, Gkik.(ti,t2) = ~iO{ti-t2){{Bi,,{h),Bl^{t2)])n, 
where B{t) arc Heisenberg operators, the braces denote 
an anticommutator, (. . .)>^ is a grand canonical average, 
and 9 is the usual step function. 

for 



Below we employ the quasiclassical method [21 



treating the inhomogeneous surface problem. In the ma- 
terials of interest a|gkj | ^ Ef for any Fermi momentum 
kj, where Ef is the Fermi energy. In addition, the su- 
perconducting energy scales (transition temperature Tc 
and the gap A) are much smaller than Ef. Under these 
conditions quasiparticles with different helicity but with 
the same k = k/|k| propagate coherently along a com- 
mon classical trajectory (determined by the Fermi sur- 
face for a = 0, ^k/ = 0), over distances much longer 
than the Fermi wavelength. We normalize gk, (gk ) = l, 
where (. . .) denotes a Fermi surface average. The quasi- 
classical propagator is then obtained as 5(k/,R, e, f) = 
f3/da/(dq)(dT)e*(i^+-)GkH.H.k-|(i + i,t-i) where 
73 is the Pauli matrix in the particle-hole space. Us- 



C/-kf/, 



and the fcrmionic anticommuta- 



mg Ly_ki^k 

tion relations for the h and h\ we derive the fundamen- 
tal symmetry relations for the 2x2 Nambu matrix com- 
ponents, 5(e,k/)22 = [(ngcr) gr(-e, -k/)ii (ngcr)]* and 
g(e,k/)2i = [(ngCr) g{-e^'kf)i2 (-icr2)(ngcr)]* {102). 



Standard procedure [21|| yields the Eilenberger equa- 
tion in helicity basis, 







(4) 



[efs - avso - A , g\+ ivfVg 

with normalization g^ = — tt^I. Here, e is the energy, 
vso = |gkf I CTafa, and A is the superconducting OP. The 



Fundamental symmetry relates 7 and 7 in helicity basis 
by7(e,k/) = [(ngcr)7(-e, -k/)(-i(T2)(ngCr)]* (1(12). In 
the bulk, 7 ~ 7+cr+ — 7!i_cr~, and 7 ~ '^"_(J^ — 7!j_<t~, 
with7^(e,k/) = -A±(k/)/(e + jV|A±(k/)|2-e2) ^nd 
7i(k/,e) = 7i(-k/,-e)*. 

The surface bound states are determined by the poles 
of the Green's function, Eq. ([6]) . We consider specular re- 
flection, whereby the component of k normal to surface 
changes sign, k ^ k, see Fig. [1] We find the amplitudes 
7k (7k), by integrating forward (backward) along incom- 
ing, k, (out goin g, k) trajectory starting from the values 
in the bulk [23|. The amplitudes 7k and 7k, in contrast, 
arc determined from the boundary conditions at the sur- 
face. Since the surface is non-magnetic, the components 
oig in the spin basis, gj. = Ul^gyfi]^, are continuous at the 
surface. This leads to a surface induced mixing of the he- 
licity bands according to C/k7kC^k = 7k = 7k = t^k7k^k' 
and Ul%U^ = 7^ = 7^ == ^k 7k^c- Fro™ Eq. dH), the 
bound states correspond to the zero eigenvalues of the 
matrix 1 — 7k 7k = 1 — 7k^k^k'7k^k^k ^* ^^"^ surface, 
and we derive our final equation for the ABS energies via 
the surface amplitudes in the helicity basis 

(l + 7+7+)(l + 7-7-) = -(l+7+7-)(l+7-7+)X • (7) 

The "mixing" factor M is determined by the change of 
gk under reflection k ^ k at the surface. 



,2 "g- 
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(8) 



where 0g, (/)g and 0g, 0g are the polar and azimuthal an- 
gles of gk and gk, respectively. If gk = gk (-B ^ B' in 
Fig. [1]) there is no helicity band mixing, M. — Q^ and we 
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FIG. 2; (Color online) Bound state energy as a function of the 
impact angle for different widths W (in units of of d//2Ao) 
of the order parameter suppression region. Blue (green) lines 
correspond to plus (minus) sign in Eq.Q. Solid lines: prin- 
cipal mode; broken lines; higher multiple reflection modes. 



recover the conditions for ABS in superconductors with 
no SO coupling. The hmit tVJ ^ oo describes pure inter- 
band scattering. In the general case {A -^ A' in Fig. [T|) 
a finite M. determines the relative weights of intraband 
and intcrband scattering. 

While assuming a uniform OP up to the surface to 
obtain the ABS spectra may seem reasonable, we show 
now that the suppression of the anisotropic (triplet) com- 
ponent of the OP in Eq. ^ near the surface drasti- 
cally modifies the ABS spectrum and the surface DOS, 
Af(e,k/) = —^ImTr {(;(£, k/)}, where Tr is a 2x2 spin 
trace, and Nf is the normal state DOS. Hereafter we 
consider a 2D material with the Rashba type SO cou- 
pling a = ajikf, gk = (k x z)/kf, and a triplet order 
parameter, A+ = — A_ = A; results for A+ ^ — A_ and 



different SO couplings will be presented elsewhere [24 1. 

To obtain insight in the role of the OP suppression, we 
consider first a simple model where A = in a layer 
of width W next to the surface, sec Fig. [2j Trajec- 
tories incident at an angle travel through a non-SC 
region of an effective width 2D = 2W/ cos (p. In this 
case Ad = cot^ 0, the surface coherence amplitudes gain 
a phase factor, 7± = 7^^ e'^eD/i^. ^ ^^ ^ ^0^ ^z2eD/vf ^ 

and the bound states are given by Im (^^e'^^'^/'"^) = 
Rc^iY+e'^^^/'"') M, which yields 



VAg 



= — tan 



2We 

Vf cos 4 



±, 



(9) 



Solutions of this equation are shown in Fig.[2l The "prin- 
cipal" modes with energies away from the continuum 
edge contribute the most to the subgap DOS. W = 
reproduces the result of Ref. |25|: each incoming trajec- 
tory yields a bound state at a different energy. For W ^ 
the main mode ei,s{4>) develops a maximum at e* < Aq, 
and we expect a peak in the surface DOS near e* due to 
abundance of trajectories contributing to N{e*). 

Fully self-consistent solution, shown in Fig. [3]^b), con- 
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FIG. 3: (Color online) Structure of the surface states for a 
Rashba triplet superconductor (A+ — — A_). (a) Order pa- 
rameter suppression; (b) DOS at locations indicated in (a); 
the large sub-gap peaks are due to the suppression of A. In- 
set: DOS for a homogeneous order parameter; note the ab- 
sence of any subgap features. (c,d) spin-resolved surface DOS 
for two trajectories; Ni{e,]i.f) (blue) and A'^x(e,k/) (green) 
correspond to blue and green branches in Fig. [2] 



firms this. Note that A 7^ at the surface, Fig. [Sj^a), as 
in other unconventional superconductors misaligned with 
respect to the interface |16|- Crucially, self-consistency 
does yield a peak in the surface DOS below the gap at 
a finite energy. Experimentally accessing this peak by 
point contact tunneling requires a sufficiently wide tun- 
neling cone as the feature arises from the trajectories at 
intermediate incident angles, see Fig. [51 

These ABS have unusual spin structure. Fig. [3][c,d) 



shows the spin resolved density of states, A'^- 



TJ, 



N ± 



N^, where A^ is the net DOS and A^"(e,k/,x) == 
— ^ImTr{cr".g(e,k/,x)}. At the interface A^^ = A^'^ = 
0. The states corresponding to different branches of 
Eq. ([9]) have opposite spin polarization. Since the spin 
polarization changes sign for reversed trajectories, the 
Andreev states carry spin current along the interface. 

Spin currents exist in NCS materials because the spin 
is not conserved, and consequently precession terms enter 
the continuity equation, atS'"(x)-hV-n"(x) = F"(x) 0]. 
Here, the spin density, ^"(x) = ^Tr/dk cr"G(k,x), the 
spin current, n"(x) == |Tr/dk {cr" , Vk}G(k,x), and 
the precession P"(x) = jrTr J dk [a"" , Vk ■ k] G(k,x), 
(where [•, •] is a commutator, and Vk = hf /m-\-afi{zx(T) 
is the band velocity), are all given in terms of Green's 
functions at imaginary relative time t = — iO. For the 
Rashba case, the precession terms are related to spin 
currents via the relations P^ = —2'mafill^, P^ = 



X 



^l 



'2maRU^, P^ = 2mafl,(H^ 
We first consider the spin currents in the normal state. 
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FIG. 4: (Color online) Spin currents near the surface. Left; 
Rashba metal. The non- vanishing components in the bulk are 



nr = -n;^ 



Near the surface Tiy is large, see text. Right: 



Spin currents and the order parameter in the SC state. 

with Ref. I3. To determine the surface spin currents we 
find the Green's function for a surface modeled as a 5- 
function barrier at a; = of strength U . The Dyson 
equation in 2x2 spin space reads G~^ = [G^'^'^]~^ — U 5{x) , 
where [G\^ ]^^ = £ ^ Ck ^ «_R(k x z)cr. For an impene- 
trable surface {U -^ oo) the solution is (for fixed A:,,)|27| 
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We solve Eq. ([TU]) numerically, and show the normal state 
surface spin currents in Fig. [H^a,b) . The most prominent 
new feature is a large surface current H^ with out of plane 
spin polarization (reminiscent to that in spin Hall bars 
[28|) that flows along the surface, and decays rapidly into 
the bulk on a Fermi wavelength scale. This component 
is related to 11^ via the continuity equation, lly{x) = 
— l/{2maii) dll^ {x)/dx. As a result, this component is 
much greater, by a factor of order kf/maji, than the bulk 
spin currents in the normal state. 

The SC spin current, shown in Fig. [Uc,d), is defined 
in the quasiclassical method relative to the normal state. 



n"-n 



N 



An;(e)(v/7V"(e,ky,x)),(ll) 



where nf(e) is the Fermi function. The surface-induced 
current with out of plane spin polarization is greater than 
the normal state current by the factor ~ TcE'i/a^ . The 
maximal amplitude at the surface is solely determined by 
the structure of the SC gap and formally survives even 
in the limit a — > 0. SC spin currents decay into the bulk 
on the scale of the coherence length, much slower than 
in the normal phase. The oscillations in Fig. [4ljc,d) are 
determined by the spin-orbit strength a and appear due 
to Faraday-like rotations of the spin coherence functions 
along quasiparticle trajectories. 



In summary, we developed a framework for the analysis 
of surface bound states and the associated spin currents 
in non-centrosymmetric superconductors, and applied it 
to a system with a Rashba-type spin-orbit coupling. We 
found that the suppression of superconductivity near the 
surface gives rise to a finite bias peak in the surface den- 
sity of states that can be probed by point contact tunnel- 
ing. We also showed that large in amplitude and slowly 
decaying spin currents with out of plane spin polariza- 
tion are carried by these surface states. This opens the 
route to future investigations of spin transport in systems 
containing superconductors without center of inversion. 
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